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Abstract 



-)— > 

'^ ' We study random networks of nonlinear resistors, which obey a general- 

Q . ized Ohm's law, V ^ F . Our renormalized field theory, which thrives on 

CJ , an interpretation of the involved Feynman Diagrams as being resistor net- 
works themselves, is presented in detail. By considering distinct values of the 

K^ ^ nonlinearity r, we calculate several fractal dimensions characterizing perco- 

[^ ' lation clusters. For the dimension associated with the red bonds we show 

^^ ■ that dred = l/i^ at least to order O (e^) , with v being the correlation length 

f^ , exponent, and e = 6 — d, where d denotes the spatial dimension. This result 



agrees with a rigorous one by Coniglio. Our result for the chemical distance, 
dj^i^ = 2 - e/6 - [937/588 + 45/49 (In 2 - 9/10 In 3)] (e/6)^ + O (e^) verifies 
a previous calculation by one of us. For the backbone dimension we find 
Db = 2 + e/21 - 172eV9261 + 2 (-74639 + 22680C (3)) eV4084101 + O (e^) , 
where C (3) = 1.202057..., in agreement to second order in e with a two-loop 



'^ \ calculation by Harris and Lubensky. 

O '• PACS numbers: 64.60.Ak, 64.60.Fr, 72.80.Ng, 05.70.Jk 
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I. INTRODUCTION 

Percolation is a leading paradigm for disorder (for a review see e.g. [|l]-3). Though it 
represents the simplest model of a disordered system, it has many applications e.g. poly- 
merization, porous and amorphous materials, thin films, spreading of epidemics etc. In 
particular the transport properties of percolation clusters have gained a vast amount of in- 
terest over the last decades. Random resistor networks (RRN) play a major role in the study 
of transport on percolation clusters for several reasons. For example, one can learn about 
the conductivity of disordered media, which might be important for technical applications. 
One can study diffusion on disordered substrates, since the diffusion constant D and the 
conductivity S of the system are related by the Einstein equation 

where e and n denote the charge and the density of the mobile particles. Nonlinear RRN, 
for which the voltage drop V over an individual resistor depends on some power r of the 
current / flowing through it, can be exploited to derive the fractal dimension of various 
substructures of percolation clusters. 

In this paper we present in detail our study of nonlinear RRN by renormalized field 
theory. A brief account of this work has been given previously in |^. It is based on an 
approach by Stephen ||^, its refinements by Harris and Lubensky [^] and its generalization 
to nonlinear resistors by Harris |0. Our work thrives on the interpretation of the involved 
Feynman Diagrams as being resistor networks themselves [§,§. This interpretation leads to a 
substantial simplification of the field theoretic calculations, as we demonstrate by calculating 
the fractal dimensions of the chemical length and the backbone to two- and three-loop order 
respectively. 

II. THE MODEL 

Consider a (i- dimensional lattice, where bonds between nearest neighboring sites are 
randomly occupied with probability p or empty with probability 1 — p. Each occupied 
bond has a finite non zero conductance a whereas unoccupied bonds have conductance zero. 
We suppose that the system is near the percolation threshold, i.e. p is close to the critical 
concentration p^ above which an infinite cluster exists. We are interested in the resistance 
Rr{x,x') between two lattice sites x and x' averaged subject to the condition, that x and x' 
are on the same cluster, 

Mr = (x(x, x')Rr{x, x'))cl{x{x, x'))c • (2.1) 

(...)c denotes the average over all configurations of the diluted lattice and xi.^^^') is an 
indicator function that takes the value one if x and x' are on the same cluster and zero 
otherwise. Note that {x{^i^'))c is nothing more than the usual correlation function in 
percolation theory. At criticality Mr obeys PJTOd 



M^ ~ |x - x'!"^"-/^ , (2.2) 

where u is the correlation length exponent defined by ,^ ~ (p — p^~" ■ 



A. KirchhofF's laws 



Here, we consider a nonlinear RRN as proposed by Kenkel and Straley |Tl]]. The bonds 
between nearest neighboring sites i and j obey a generahzed Ohm's Law, 

(2.3) 





j i Phj ij ij\ ' 


or equivalently. 






^^AVJ-V^\VJ-V^r'=kJ, 



(2.4) 

where aij (pij) is the nonhnear conductance (resistance) of the bond {i,j), lij is the current 
flowing through the bond from j to i and Vt is the potential at site i. The exponents r and 
s are describing the non linearity with r = s~^. The conductance and the resistance are 
related via aij = p~j. 

Suppose a current / is put into a cluster at site x and taken out at site x'. Those sites 
connected to x and x' by mutually avoiding paths are constituting the backbone between x 
and x'. The power dissipated on the backbone is by definition 

p = /(v;-v;o . (2.5) 

Using Ohm's law, it may be expressed entirely in terms of voltages as 

P = Rr{x,x')-' IK - V;.r+' = Yl ^^M - Kf ""' = P{{V}) , (2.6) 

{id) 

where the sum is taken over all nearest neighbor pairs on the cluster and {V} denotes the 
corresponding set of voltages. As a consequence of the variation principle 



_d_ 



\pi{v})-T.W^ 



' + - 



, (2.7) 



one obtains the circuit equations 



y: ct,, m - V,) \v, - v,r' = -y: u, = k , (2.8) 

0-> {3) 

where li = I {Si^ — (^i.^O and the summations extend over the nearest neighbors of i. 
Alternatively to Eq. ([^.6|) the power can by rewritten in terms of the currents as 



P = Rr{x,x'm■+' = YP^M^J^' = PiU}) , (2-9) 

{id) 

with {/} denoting the set of currents flowing through the individual bonds. Obviously 
the cluster may contain closed loops as subnetworks. Suppose there are currents I /^'^ \ 
circulating independently around a complete set of independent closed loops. Then the 
power is not only a function of I but also of the set of loop currents. The potential drop 
around closed loops is zero. This gives rise to the variation principle 

5^P ({/'"}. /)=0. (2.10) 



Eq. (|2.10| ) may be used to eliminate the loop currents and thus provides us with a method 



to determine the total resistance of the backbone via Eq. (|2.9| 



B. Connection to cluster properties 

Here we provide background on the meaning of (pr for some specific values of r. For 
r — > 1, one recovers the hnear RRN. 0i is the usual resistance exponent as studied to order 
e\ e.g., in §. 

Other values of r are related to the fractal dimension of substructures of percolation 
clusters. Consider r -^ — l"*". One obtains immediately as a consequence of Eq. (|2.9| ), that 



i?_i(a;,x')= lim Vpij 



r-+-l+ 



(i,i> 



J r+1 



J2 Pi,J ' 



(2.11) 



with only those bonds carrying non zero current contributing to the sum on the right hand 
side. Hence 



M_i(x,x') 



Mb 



(2.12) 



where Mb stands for the number of bonds belonging to the backbone. Thus, the fractal 
dimension of the backbone can be expressed as 



Db = lim (pr/^ 



(2.13) 



Now we turn to r ^ oo and r — > 0^ following the lines of Blumenfeld and Aharony 
T^ . On the backbone between two sites x and x' one may distinguish between two different 



substructures: blobs formed by multi-connected bonds and singly connected bonds which 
are referred to as red bonds. Both substructures are contributing to the resistance of the 
backbone 



blob 
Rr{x,x') = J2 Pid 
(id) 



T r+1 red 

+ 1^ P^,] ' 



(2.14) 



where the sums are taken over all bonds belonging to blobs and over all red bonds respec- 
tively. Since sites on a blob are multi-connected by definition, /j ., < /, and thus 



blob 



lim ^ pi^j 



{i,i) 



^i,j 



r+1 



. 



(2.15) 



In conclusion, the dimension of the red bonds is related to d>r via 



dj-ed = lim (/)j./z/ . 



(2.16) 



Consider now the first site x at some end of a blob. An entering current / splits into 
currents li^^ flowing to nearest neighbors i with 



0"j,x IK - ^j| 



(2.17) 



In the limit s ^ oo the ratios |/j,x|/|-^j,z| vanish whenever ai^^ |Kc ~ ^T < ^j,x \Vx — Vj\^- 
Thus, current flows only through the resistor with the largest Uj^^ |\4 ~ ^T- This argument 



may be iterated through the entire blob. One identifies either a single self avoiding chain 
through which I fiows, with 



{i,3) 



(2.18) 



being the power dissipated on the chain, or several of such chains with identical power. 
The expression in Eq. ( |2.18| ) is minimal for minimal J2(i,j)Pi,jy i-e., the current chooses the 
shortest path through the blob and one is led to 



dmin = lim 0,./z/ 



(2.19) 



for the chemical length exponent. 



C. Generating function 

Our aim is to determine M^. Hence our task is twofold: we need to solve the set of 
Kirchhoff 's equations ( p.8|) and to perform the average over all configurations of the diluted 
lattice. It can be accomplished by employing the replica technique [§],^. The network is 
replicated D-fold: I4 ^ 14 = (v^^\ . . . , V^^y One introduces 

tPxix) = exp (tX ■ V,) , (2.20) 

where X -Vx = Y.a X^"^V^"^ and A 7^ 0. One considers the correlation function 

G {x,x'; X) = (^x(a^)V^-A(^'))_ (2.21) 

given by 

1 



rep 



G (x, x'; A) = /z-^ I Yin dVt^ exp [--l^P {{v}) + y E ^/ + ^^ ■ {t - V. 



j 0=1 



. (2.22) 



c 



Here P({f}) = ELi P {{v^''^}) = Ef=i E(,,) a.,, 






s+1 



ization 



z=/nrf^,exp[--i^p(m) + ^EK^ 



and Z is the normal- 



(2.23) 



Note that we have introduced an additional power term ^ X^i ^i'- This is necessary to 



give the integrals in Eqs. ( |2.22| ) and ( |2.23| ) a well defined meaning. Without this term the 
integrands depend only on voltage differences and the integrals are divergent. Physically 
the new term corresponds to grounding each lattice site by a capacitor of unit capacity. The 
original situation may be restored by taking the limit of vanishing frequency, a; — > 0. 

In contrast to the linear network, P is not quadratic and hence the integration over 
the voltages is not Gaussian. This obstacle may be surmounted by employing the saddle 
point method 0. The saddle point equation is identical to the variation principle stated in 



Eq. (|2.7|) . Thus the inaxiinum of the integrand is determined by the solution of the circuit 
equations (p.8|) and, up to an unimportant multiphcative constant which goes to one in the 
hmit D -^0, 

G {x, x'; a) = /exp [ -^^Rr (x, x') J \ (2.24) 

where 

A. (a) = E (-.A("))^''' . (2.25) 



a=l 



Consequently, G f x, x'; Aj may serve as a generating function for M^, which may be obtained 
by taking the derivative of 

G[x,x';X) = {xix,x'))cU + ^^Mrix,x') + ...] , (2.26) 

with respect to A,, evaluated at A^ = 0. 

At this point a comment on the nature of A is appropriate. One sets 

X(») = zAo + e^°) , (2.27) 

with real positive Aq and ^^°''> and imposes the condition Y.a=i C''"^ = 0. The saddle point 
approximation in Eq. ( ^.24| ) may be justified by demanding 

Ao > 1. (2.28) 

On the other hand, substitution of Eq. (|2.27|) into the definition of A^ leads to 

A. (A) = E |aS+^ -^{r + l) A^^^") - !±±}lx^,-^^(»)^ + . . .| 

= DXl-^' - '-^^y^^e + ... . (2.29) 

Thus one can justify the expansion in Eq. ( ^.26[ ) by invoking the conditions 

AS+^ < D-^ and X^-^^ < 1 . (2.30) 

Note that the replica limit D ^ allows for a simultaneous fulfilment of the conditions 
( p.28| ) and ( p.30| ). However, we will not only rely on these conditions on A. We will provide 
several consistency checks for the validity of Harris' saddle point approach as we go along 
and reproduce known results. 



D. Field theoretic Hamiltonian 

Since infinite voltage drops between different clusters may occur, it is not guaranteed 
that Z stays finite, i.e., the limit lim/j^o ^^ is not well defined. Moreover, A = has to be 
excluded properly. Both problems can be handled by restoring to a lattice regularization of 
the integrals in Eqs. ( p.22|) and ( p.23| ). One switches to voltage variables Q = A6k taking 
discrete values on a D-dimensional torus, i.e. k is chosen to be an D-dimensional integer 
with -M < A;(°) < M and fc(°) = fc(°)mod(2M). In this discrete picture there are (2M)^-1 
independent state variables per lattice site and one can introduce the Potts spins ||13[ 



$. 



{X 



(2M) 



-D 



Y^ exp {iX ■ e) i)^{x) 



h.e. - (2M) 



-D 



(2.31) 



subject to the condition Y,e^e{x) = 0. 

Now we revisit Eq. ( |2.22| ) . Carrying out the average over the diluted lattice configurations 
provides us with the weight exp(— if^ep) of the average (■ ■ ■)rep, 



H, 



rep--ln(exp(^--^P(^-) + ^E^^ 



0^ 



(iJ) 



J2 In [exp [-—^(^i,j \di - dj 



s+l 



C 



C 



lUJ 



E« 



i ) 



(2.32) 



where we have introduced the abbreviation 16' 



s+l 






^(") 



s+l 



By dropping a constant 



term Nb ln(l —p), with Nb being the number of bonds in the undiluted lattice, one obtains 



{hi) i 



(«.i> el' * e 



(2.33) 



where 



K{e) 



InU 



P 



1 — p 



exp 



s + l 



-a 



|S + 1 



(2.34) 



and 



'=(»>fE»7 



(2.35) 



Note that K (9j is an exponentially decreasing function in replica space with a decay rate 
proportional to a^^. For large a, the Hamiltonian Hj-^p describes a translationally invari- 
ant short range interaction of Potts spins in real and replica space with an external one 
site potential hiOj. Moreover, the interaction potential K (9j is an analytic function of 






OM 



s+l 



Thus the Fourier transform 



k (a) = ^ X: exp (-.A ■ e) K {e) (2.36) 

can be Taylor expanded as 

oo 

k(\)=WQ-Y. wr,p [-A. (a)]" , (2.37) 

with Wo and Wr,p ~ cr"^ being expansion coefficients. 

In the hmit of perfect transport, a -^ oo, K \9\ goes to its local limit K \9\ = K5g^, 
with K being a constant. The interaction part of the Hamiltonian reduces to 

HZ = -KT.T.^ei^)^ei3) ■ (2-38) 

This represents nothing more than the (2M) states Potts model which is invariant against 
all (2M) ! permutations of the Potts spins $g-. If a^^ ^ 0, this S.^m)^ symmetry is lost in 
favor of the short range interaction. 

We proceed with the usual coarse graining step and replace the Potts spins $g'(x) by 
order parameter fields ip ( x, 9) which inherit the constraint J^ff^ (x, ^ = 0. We model 
the corresponding field theoretic Hamiltonian Ti in the spirit of Landau as a mesoscopic 
free energy from local monomials of the order parameter field and its gradients in real and 
replica space. The gradient expansion is justified since the interaction is short ranged in both 
spaces. Purely local terms in replica space have to respect the full S,2M)° Potts symmetry. 
After these remarks we write down the Landau-Ginzburg- Wilson type Hamiltonian 

r/-^2 Wr / ::\ J^ f d 






dei") 




+ ^(Vv.(x,^-)) +^^(x,^-) +^^-^^(x,^-)|. (2.39) 

Here we have neglected all terms that are irrelevant in the renormalization group sense, r 
and Wr are now coarse grained analogues of the original coefficients wo and Wr,i appearing 
in Eq. ( 2.371 ). Terms associated with Wr^p are irrelevant for p > 2 and therefore neglected. 



Note again that 7i reduces to the usual (2M)^ states Potts model Hamiltonian by setting 
w^ = as one retrieves purely geometrical percolation in the limit of vanishing Wr (cr -^ oo). 

III. RENORMALIZATION GROUP ANALYSES 

A. Resistance of Feynman Diagrams 

The diagrammatic elements contributing to our renormalization group improved pertur- 
bation calculation are the three point vertex —g and the propagator 



'^A.O 1 '^A,0 



p2 + r-«;,.A^(A) p2 + r - w^A,. (A) P^ + r 



(3.1) 



Note that we have switched to a (p, A j -representation by employing Fourier transformation 
in real and replica space. Eq. ( |3.1| ) shows that the principal propagator decomposes into a 
propagator carrying A's (conducting) and one not carrying A's (insulating). This allows for 
a schematic decomposition of principal diagrams into sums of diagrams consisting of con- 
ducting and insulating propagators (see App. 0). Here a new interpretation of the Feynman 
diagrams emerges [Q. They may be viewed as resistor networks themselves with conduct- 
ing propagators corresponding to conductors and insulating propagators corresponding to 
open bonds. The parameters s appearing in a Schwinger parametrization of the conducting 
propagators, 



1 



p^ + r — WjAr ( A 



ds exp —s (t + p^ — WrAr (a 



(3.2) 



correspond to resistances and the replica variables iX to currents. The replica currents are 
conserved in each vertex and we may write for each edge i of a diagram, Aj = Aj (A, {/?}), 

where A is an external current and {k} denotes a complete set of independent loop currents. 
The A-dependent part of a diagram can be expressed in terms of its power P: 



exp Wr ^ SiAr (Xij j = exp WrP (a, {k} 



(3.3) 



The new interpretation suggests an alternative way of computing the Feynman diagrams. 
To evaluate sums over independent loop currents 

(3.4) 



^exp \wrP [X,{k,} 

{k} 



we employ the saddle point method under the conditions discussed at the end of Sec. pC 



Note that the saddle point equation is nothing more than the variation principle stated in 
Eq. ( |2.1CI| ). Thus solving the saddle point equations is equivalent to determining the total 
resistance R{{si}) of a diagram, and the saddle point evaluation of ( p.4|) yields 



exp 



Rr {{Si}) WrAr (A 



(3.5) 



where we have omitted once more multiplicative factors which go to one for D ^ 0. A 
completion of squares in the momenta renders the momentum integrations straightforward. 
Equally well we can use the saddle point method which is exact here since the momentum 
dependence is purely quadratic. After an expansion for small A^ (Xj all diagrammatic 
contributions are of the form 



/ (p2, A^) = Ip (p2) + I^ (p2) ^^A^ (a) + . . . 

n dSi [l + Rr ({S,}) WrAr (x) + . . .] D {p\ {sj 



(3.6) 



D (p^, {si}) is nothing more than the integrand one obtains upon Schwinger parameteriza- 
tion of the corresponding diagram in the usual (p^ theory. 



B. Renormalization and scaling 



We proceed with standard techniques of renormahzed field theory [|I?|. The ultraviolet 
divergences occurring in the diagrams can be regularized by dimensional regularization. We 
employ the renormalization scheme 



^ -. ^ = z^/2^ 



f = z-^z^ 



Z Zyj^Wj. , 



g^g = Z-'/'Zi/'G:'/\''-^^ 



/2,,^/2 



(3.7a) 
(3.7b) 



where e = 6 — d and /i is an inverse length scale. The factor G^ = (47r)^'^/^r(l + e/2), 
with r denoting the Gamma function, is introduced for convenience. The Z factors may be 
determined by minimal subtraction, i.e., they are chosen to solely cancel poles in e. Z, Z^ 
and Zu are the usual Potts model Z factors. They have been computed to three loop order 



by be Alcantara Bonfim et al |25|. It remains to calculate Z^^. We postpone this calculation 
to Sec. |IV|. 

The unrenormalized theory has to be independent of the length scale ^~^ introduced 
by renormalization. In particular, the connected A^ point correlation functions must be 
independent of /i, i.e.. 



d 



f^g-GN ({p, Wr Ar (a) } ; T, ^) =0 

for all A^. Eq. (|3.8|) translates via the Wilson functions 



(3.8) 



p{u) = fl 



du 



Criu) 



dlnwr 



k{u) = jl 



d\n.T 



fx- 



dfi 



d 



, 7... (m) =/i— InZ. 



(3.9a) 
(3.9b) 



where the bare quantities are kept fix while taking the derivatives, into the Gell-Mann-Low 
renormalization group equation 



' d ^d d a N ' 



du du 



dr 



dWr 



G'jv({x,w,.A,. (A)};r,M,/i) =0 . (3.10) 



The particular form of the Wilson functions can be extracted from the renormalization 
scheme and the Z factors. 

The renormalization group equation is solved by the method of characteristics. At the 
infrared stable fixed point u*, determined by /3 (u*) = 0, the solution reads 



G'jv({x,u7,A, (A)};r,u,/i) = P*^/^Gn {{l^J^'wrAr {x)} ■J'''T,u*,lfi) , 



(3.11) 



where 7* = 7 (m*), k* = k {u*) and (* = (r (u*). 

To get a scaling relation for the correlation functions, a dimensional analysis remains to 
be performed. It yields 

Gn ({x, WrAr (a) };r,u,fi)= /i('^-2)^/2G'^ ({/ix, fi-'^Wr-Ar (X) } ; /i-V, u, 1) . (3.12) 



10 



From Eqs. (|3.11|) and ( p.l2| ) we drive the scaling relation 

Gn ({x, WrK (a) } ; r, m, /i) = /('^-2+r?)JV/2^^ (^^i^^ l-'I'^/'^WrK (a) } ; l-^'^'r, u*, n 
with the well known critical exponents for percolation [O 



1 



206 



V = l 



21 9261 



e^ + 



93619 256 



8168202 7203 



C(3) 



e' + OU 



and 



z/ = (2 - K*) 



_i _ 1 5 589 2 

~ 2 ^84^^ 37044^ ^ 



716519 



89 



130691232 7203 



C(3) 



.^ + O ie' 



(3.13) 



(3.14) 



(3.15) 



Note that C, in Eq. ( p.l4| ) stands for the Riemann zeta function and should not be confused 
with the Wilson function defined above. The exponent 0^ is defined by 



(t)r = v{2-C) = ^{'^-V + ^r) 



(3.16) 



with ipr = 7«,r ('^*)- It has been calculated for arbitrary r to one-loop order by Harris [0. 
He found 



1 + -^c. + O (e^) 



where 



1 /■! 



d^ 



:i-a 



1 + 0'/'^ + (1-0 



l/r 



(3.17) 



(3.18) 



Our calculation gives the same result. 

Equation ( ^.13| ) implies the following scaling behavior of the two point correlation func- 
tion G = 6*2 at criticality, 



G ( |x - x' \;wr Ar ( A) ) = I'^-^+'^G (/|x - x'l; /-'^'■/^w^A, (A 



(3.19) 



where we dropped several arguments for notational simplicity. The choice / = |x — x'| ^ and 
a Taylor expansion of the right hand side of Eq. ( p.l9| ) lead to 

G (|x - ^'\;WrAr (a)) = |x - xf -'^-'' (l + w,A^ (a) |x - x'l'^'-/^ + . . .) . (3.20) 

Comparison with Eq. ( p.26| ) gives us the scaling behavior of the average resistance: 



MJx,x')~ Ix-x'l*'-/^ 



(3.21) 



IV. FRACTAL DIMENSIONS 



In this section we calculate (hr for r — * oo, r ^ 0+ and r 



— 1+. As discussed in 

Sec. [II B| , this provides us with the fractal dimension of the red bonds, the chemical length, 
and the backbone respectively. 
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A. Red bonds 

Consider r -^ oo. As argued in Sec. [11 B|, blobs do not contribute to the total resistance. 



Now we take direct advantage of our view of the Feynman diagrams as being resistor networks 
themselves. In analogy to real networks, the resistance of closed loops vanishes. Only singly 
connected conducting propagators contribute to the total resistance of a diagram, i.e. 



singly 

E 



^oo i{s^}) = E ^* ' (4-1) 



with the sum being taken only over singly connected conducting propagators. The contri- 
bution of a diagram to the renormalization factor Z^^ takes the form 

, , roc singly 

iw (p^) = / n^^. E '^D (p^ {s,}) . (4.2) 

° J i 

Note that a factor Sj in Eq. ( [4.2|) corresponds to the insertion of |y9^ into the ith. edge of 
the diagram. We generate Iw (p^) by inserting ^{p"^ in each singly connected conducting 
propagator. This procedure is carried out up to three loop order, i.e., every conducting 
propagator in App. ^ that does not belong to a closed loop gets an insertion. For detail see 
App. 0. The resulting diagrams are displayed in Fig. |l|. 

Now consider the contributions of the diagrams listed in App. |^ to Zr- These can be 
generated by inserting ^ip"^ in conducting as well as in insulating propagators. Again, one 
obtains the diagrams depicted in Fig. |l| with the same fore-factors. Consequently, Z^,^ and 
Zr are identical at least up to three-loop order. The same goes for the corresponding Wilson 
functions Coo and k. From the definition of (pr it follows that 

2 — t* / 



2- K* 



1 + e^ . (4.3) 



Note that this result is in agreement with the rigorous one by Coniglio |16,17|, who proved 
that (ired = l/'^- We rate this as an indication for the validity of Harris' saddle point 
approach. 

B. Chemical length 

In the limit r — i> 0"*" only the shortest self avoiding path of conducting propagators 
contributes to the total resistance of a diagram. In other words, the total resistance has to 
be determined such that 

T. T. ^^ (4-4) 

paths igpatli 

is minimal, where the first sum is taken over all self avoiding paths of conducting propagators 
connecting the external legs of a diagram. 

For r ^ 0+ the conducting propagator reads 
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r + p2 + iwo Ea=i ^^"^ X + iwoX ' 



(4.5) 



with X = T + P^ and A = I]«=i '^''"^- We start with the two one-loop diagrams A and B (see, 
App. ^. The diagram A translates into 

A = — - / / dsids2 exp [siXi + S2X2] exp [-iwoX min (si, S2)] 

Z Jp JO 

/• /"OO 

= g^ I I dsids2 exp [siXi + S2X2] exp [-iwoAsi] 9 (s2 - si) , (4.6) 

Jp JO 

where 9 denotes the step function and /p is an abbreviation for (27r) / d^p. Diagram B reads 

^ = IT / / dsids2e^V [^iXi + 52X2] exp [-iwqXsi] , (4.7) 

2 Jp JO 

and hence, 

k-2B = -g^ I I (isirfs2exp [siXi + S2X2]exp [-■ iwoAsi] 6'(si - S2) . (4.8) 

Jp JO 

Now we take a short detour and present some features of the field theory of dynamical 



percolation as studied by one of us some time ago |TB|. The dynamical functional J that 
leads to the diagrammatic expansion for the calculation of correlation and response functions 
is given by 



J= f d'^xdt-fip 



1 



d 



-^^ + (--A)-.w-2 



9 ~ 



^ 



Here, (f) (x, t) = 7 J_^ dt'if (x, t') and (f (x, t) is the response field. 

Gi,i(x,t) = (^(x,t)(^(0,t))^ 



(4.9) 



(4.10) 



is the density response function that describes a growing cluster initiated by a germ at 
(x = 0, t = 0) which percolates at the critical point. Near this percolation point the response 
function scales as 



G'i,i(x,t) = e 



-{d-2+ij) 



f (x/e, t/e 



(4.11) 



where f is a scaling function, ^ = \t\ '^ is the correlation length, and z is the dynamic 
exponent given to second order in e in [^ . 



The diagrammatic elements of dynamical percolation are the propagator 

G (p, t)=9 (t) exp [-7 (r + p^) t 



(4.12) 



and the vertices 7(7 and —'yg9 {t — t') . These elements are depicted in Fig. 0. Note that the 
Fourier transformed propagator reads 



G{p,u;) 



1 

icij + 7 (r + p2 



(4.13) 
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and can be identified with ([4.5|) up to a factor 7"^ by setting •ywoX = u. Thus, the renor- 
mahzation of Wq is directly related to that of the kinetic coefficient 7. One finds that z is 
related to the chemical length dimension by 2; = (imin- 

The one loop contribution to the vertex function Fii (p,u;) = l/Gi 1 (p, cj) is visualized 
in Fig. |. We find 

a = - {igf I I dtidt2 exp [7 (tiXi + ^2X2)] (ti - t2) exp [-iujti] 
Jp Jo 

= A-2B (4.14) 

if we identify 'yU = Sj. 

Now we turn to the two-loop diagrams. In the same manner as in the one-loop case we 
obtain 

C - 4D - E + 2F + 4G 

moo J' / ^ \ 

n dsi exp J2 ^iXi 

X \ exp [-IwqX (si + S4 + S5)] \6 (s2 - si - S5) 61 (53 - S4 - S5) 

- 6 (s2 - si - S5) - 6* (s3 - S4 - S5) + 1} 

+ exp [-iwqX (si + S3)] \9 (s2 -h S4 - si - S3) 9 (s4 + S5 - S3) 9 (s2 + S5 - si) 

- 6* (S2 + S4 - Si - S3) - 61 (S4 + S5 - S3) - 6^ (S2 + S5 - Si) + 2 U 

= (7^)^ / / / X{dtie^^[Y^ tilXi 

JpJciJO ^^^ \i^^ ) 

X I exp \-i^ {ti +U + ts)] 6 {h + ts - h) 9 {U + ts - ^s) 

+ exp [~iuj (ti + t3)] {d (ii - ^2 - ts) ih + h- U) 

+ 9{h+t3-t2- U) [9 {U -h-h) + 9 (t3 - t4 - h)] }} 

= b + c + d + e. (4.15) 

The diagrams b, c, d, and e are depicted in Fig. ^. For the second bold two-loop diagram 
we find 

H - I - 2J + 2K + L 

moo 5 / 5 \ 

n dsi exp Yl ^iXi 

X I exp [-iwoX (si + S2 + S3)] \^9 (S5 - si - S2 - S3) 9 {34 - S2) 

- 9 (S5 - si - S2 - S3) - 9 (S4 - S2) + 1} 

+ exp [-iwoXs5] -{9 (si + S2 + S3 - S5) 9 (si + S4 + S3 - Sg) 

- 6^ (si + S2 + S3 - S5) - 6^ (si + S4 + S3 - S5) + l| 
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= (7^)^ / / / HdUexpiY^ UlXi 

X I exp [-iu (ti + t2 + ta)] (ta - ^4) (ti + ta + ^3 - ^5) 

+ exp [-itoh] 9 (ts - ti - t2 - tg) ^ (t2 - t4) } 

= f + g. (4.16) 

The diagrams f and g are shown in Fig. ^ 

The dynamic diagrams lead to the result for the dynamic exponent z stated in |]18 . 
Since we identified the two diagrammatic expansions up to two-loop order, the RRN gives 
the same result for the chemical length dimension as the dynamic approach in | |18|| , 



d =2--- 

O 



937 45 / „ 9 . \1 /e^2 



588 49 V 10 



— In 2 In 3 



- +0(6^) . (4.17) 



Moreover, another consistency check for the saddle point approximation is fulfilled. 

Obviously, dmin has to approach one for d ^ 1. This feature can be incorporated by a 
rational approximation yielding 

Cimin ^ 1 + (1 - (1 + 1^ - 0.03016^) . (4.18) 

Due to the rich structure of rj in the percolation problem, 

3 365 + 1401n2-1261n3 2 ...^^ 

^0 = e e (4.19) 

^" 14 5488 ^ ^ 

might be better suited for such a comparison than t/min. It is known exactly that ipo vanishes 
in one dimension. A rational approximation yields 

V'o ~ (1 - (-^e - 0.1018e2) . (4.20) 



drain and ipo are compared to numerical simulations by Grassberger |]T9[ in Fig. ||. The 



rational approximants agree reasonably well with the numerical estimates at d = 3. At 
d = 2, the approximant for dmin seems to be in conformity with the simulation result. 
However, the good agreement should be taken with caution. It might be accidental, since 
-00 {d = 2) hardly agrees with the numerical value. 

C. Backbone 



Now we focus on the limit r -^ — 1"*". As argued in Sec. PB| , the resistance of the 
backbone between two sites x and x' is given by 

i?_i(x,x') = EPM' (4-21) 

with the sum running over all current carrying bonds of the underlying cluster. In analogy, 
the resistance of a Feynman diagram is given by 
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cond 



(4.22) 



where the sum is extending over all conducting propagators of the diagram. The contribution 
of a diagram to Z^_j now takes the form 



/•oo cond 

^ ■ n 



(4.23) 



We proceed in the same manner as in Sec. [IV A|. However, now ^0^ is inserted into all 



conducting propagators. For details of the calculation see App. |C[ Minimal subtraction 
leads to the renormalization factor 



^UI_l 



i + — + 



4e 



12 144 



'^-|C(3). 



Oiu 



Via the Wilson functions we obtain the exponents 

2 r 2075 



^- 



1 



' ' 7 • + 



16C(3)- 



126 



i +o(^') 



and 



D 



B 



2 + — e 
21 



172 



-e^ + 2- 



-74639 + 22680C (3) 



+ 0(t 



(4.24) 



(4.25) 



(4.26) 



9261 4084101 

Note that our result agrees to second order in e with calculations by Harris and Lubensky 



20| based on another approach. This is again in favor of the saddle point approximation. 
In Fig. 1^ we compare the e-expansions as well as the rational approximants 



ij. 



2e^ f t\ f t 

1- - 1 + 1.2625 

49 V 5/ V 500 



and 



D 



B 



5 



26 
105 



e + 



7166 
231525' 



0.0170e^ 



(4.27) 



(4.28) 



to numerical simulations by Grassberger ||2T| and Moukarzel p2[. At d = A the results 
agree within the numerical errors. However, a higher accuracy of the numerical estimate is 
desirable. At d = 3 and d = 2 the analytic results look less realistic, but they reproduce the 
shape of the dependence on dimensionality. 



V. CONCLUSIONS AND OUTLOOK 

By employing a saddle point approach due to Harris we calculated the exponent (pr/i^ for 
the critical behavior of the resistance in a diluted network. We focussed on distinct values 
of the nonlinearity r, namely those related to the fractal dimensions of the red bonds, the 
chemical path and the backbone respectively. 
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We provided several consistency checks for the saddle point approach. The validity of 
the approach seems to be beyond question. 

For dimensions close to the upper critical dimension six, our results for rfmin and Db 
are the most accurate analytical estimates that we know of. The analytic results agree 
reasonably well with the available numerical simulations. At low dimensions the agreement 
becomes less pronounced. 

Our interpretation of Feynman diagrams proved to be a powerful tool. It simplified the 
renormalization group improved perturbation calculation considerably. The technique used 
here may be applied to other aspects of transport on percolation clusters. For example it can 
be employed to calculate the family of noise exponents {ipi} for diluted resistor networks, 
as treated by Park, Harris and Lubensky |^ to one-loop order. Our two-loop calculation 
yielding 



e^ 



"^^ ^^7(1 + 0(1 + 2/) ^ 12348 (1 + if (1 + 2/)' 

X I 313 - 6727 + /|3327 - 40327 - 8/{4 (-389 + 2737) 

+ / [-2076 + IOO87 + / (-881 + 3367)] }| - 672 (1 + if (1 + 2lf ^(1 + 21) I (5.1) 

will be reported in a separate publication in the near future. In Eq. ( p.l| ) 7 denotes Euler's 
constant and \1/ stands for the Digamma function. 
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APPENDIX A: DECOMPOSITION OF DIAGRAMS 

Here we list the decomposition of the primary two leg diagrams (bold) into conducting 
diagrams composed of conducting (light) and insulating (dashed) propagators. The listing 
extends up to three-loop order. Note that the conducting diagrams inherit their combinato- 
rial factor from their bold diagram. For example, the diagrams A and B introduced below 
have to be calculated with the same combinatorial factor, namely \. 
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APPENDIX B: DIAGRAMS FOR THE RED BONDS 

This appendix gives details on the diagrammatic contributions to the renormahzation of 
Wr and r in the hmit r ^ oo. As an example we consider the one-loop diagrams A and B. 
In Sec. |IVA| we argued, that only singly connected conducting propagators contribute to 



Thus, A gives no such contribution at all. The contribution of B can be expressed as 



(Bl) 



where the lines stand for conducting propagators evaluated at zero currents and the solid 
dot for an ^0^-insertion. The resulting contribution of A — 2B is 




-2 




(B2) 



Now we turn to r. Z^ can be calculated by inserting ^0^ into conducting and insulating 
propagators. The contribution of both A and B reads 




(B3) 



The resulting contribution of A — 2B is again the one stated in ( P2|) . 

We carry out the insertion procedure for both, Woo and r, up to three-loop order. One 
obtains in both cases the same diagrams with the same fore-factors. The result is listed in 
Fig. 1. 



APPENDIX C: EVALUATION OF DIAGRAMS FOR THE BACKBONE 

In this appendix we give some details of the calculation of the backbone dimension. As 
described in Sec. p.V C] , we insert |0^ into each conducting propagator. The diagram A for 
example has two conducting propagators. Its contribution to Z^_^ can be expressed as 
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(CI) 



As in App. P, the lines stand for conducting propagators evaluated at zero currents and the 
solid dot for a i0^-insertion. The diagram B contributes via 




(C2) 



and hence the total contribution of A — 2B vanishes. 

The procedure is carried out up to three-loop order. It results in 







(C3) 



These diagrams can be most conveniently evaluated by mapping them onto those calculated 



in [|T5|. The two-loop contribution for example can be re-expressed as 





(C4) 



Note that we have explicitly extracted the combinatorial factor ^ from the diagram. This is 
important at this stage, since different diagrams, each inheriting its combinatorial factor from 
its bold diagram, may be mapped onto the same three-leg diagram. The additional factor 
on the right hand side cancels the combinatorial factor | of the three-leg diagram as well as 
a vertex —g. Similar identifications can be made for the three- loop diagrams appearing in 
Eq. (p3|) . After all, the following diagrammatic contributions to the renormalization of W-i 
are obtained: 



+ : 





(C5) 
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where we have dropped an overall factor ~g~^. 

The e-expansion results for the diagrams in Eq. (|C5| ) can be gathered from [15]. However, 
we did not entirely rely on the results stated there. We also did the calculations on our own 
and found the same results leading to the renormalization factor given in Eq. ( |4.24| ). 
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FIGURES 



FIG. |T] 

Diagrammatic expansion in the limit r — > oo. The listed diagrams including their fore-factors 
can be obtained from the conducting diagrams shown in App. in two different ways: firstly, 
by inserting ^ip"^ into all singly connected conducting propagators and secondly, by inserting 
|y9^ into every conducting and insulating propagator. As a consequence, the renormalization 
factors Zw^ and Z^ are identical. The lines stand for conducting propagators evaluated at 
zero currents, the solid dots for iy^^-insertions. 



FIG. I 

The propagator G (p, t) as well as the vertices 'jg and —'jgO [t — t') (from left to right). 



FIG. I 

In the limit r -^ 0^ we map the bold one-loop diagram (see App. ^ onto the dynamic one 

shown here. The meaning of the graphic elements may be inferred from Fig. 0. 



FIG. I 

Dynamic diagrams obtained in the limit r 



FIG. I 

Dynamic diagrams obtained in the limit r 



FIG. I 

Dependence of the exponents dmm and ipo = dmin — l/^ on dimensionality. The e-expansion 
(full squares) and the rational approximation (open squares) are compared to numerical 
simulations (circles). For d^i^i we take Monte Carlo results by Grassberger. At d = 2 we 



insert the exact values @,3 i' = 4/3 and 7 = 43/18. At rf = 3 we use Monte Carlo results 



by Ziff and Stell PSI: z/ = 0.875 ± 0.008, 7 = 1.795 ± 0.005. 



FIG. 

Dependence of the exponents Db and ip^i = Db — ^jv on dimensionality. The e-expansion 
(full squares) and the rational approximation (open squares) are compared to numerical 
results (circles) by Grassberger id = 2) and Moukarzel [d = 3,4). They determined Db by 
simulations. At d = 2 and d = 3 we use the same values for 7 and z/ as in Fig. ^ At d = 4 



we take z/-^ = 1.44 ± 0.05 [||] and 7 = 1.44 | 



25 




+2 



-4 



-4 



-20 



-4 



-2 



-3 



+ 














-20 



-2 








+2 



-4 








-16 






-2 







-2 



-4 



-2 













FIG. 1. 



26 



f) < CO 




t / 



FIG. 2. 




FIG. 3. 




+ 





+ 




FIG. 4. 




+ 




FIG. 5. 



27 




FIG. 6. 



28 



1 


2 


3 


4 5 


2.01 




- s 




_ 








^4\ 


& 
-^^ 


o- 


1.5 






^ 




1.0 








1 


0.5 

n n 








\ 

\ 



2.0 



1.5 



1.0 



0.5 




FIG. 7. 



29 



